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' We use analytic perturbation theory to present a new approximate metric for a 

rigidly rotating perfect fluid source with equation of state (EOS) e + (1 — n)p = eo- 
Q I This EOS includes the interesting cases of strange matter and Wahlquist's metric. It 

Q I is fully matched to its approximate asymptotically flat exterior using Lichnerowicz 

■ junction conditions and it is shown to be a totally general matching using Darmois- 

Israel conditions and properties of the harmonic coordinates. Then we analyse the 
Petrov type of the interior metric and show, first, that in the case corresponding to 
Wahlquist's metric it can not be matched to the asymptotically flat exterior. Next, 
that this kind of interior can only be of Petrov types I, D or (in the static case) O 



bi)' and also that the non-static constant density case can only be of type I. Finally, we 

check that it can not be a source of Kerr's metric. 

> '. PACS numbers: 04.25.Nx, 04.40.Dg 

in 

^ ■ I. Introduction 

(N 

One of the regrettable facts about General Relavity is that, by now, despite 
the numerous exact solutions and modern methods to generate them, it has failed 
^ ■ to provide an exact solution describing a rotating stellar model, i.e., a spacetime 

^ . corresponding to an isolated self-gravitating rotating fluid in equilibrium. These 

stellar models are built matching an interior spacetime describing the source and the 
exterior spacetime that encloses it. A candidate interior solution should correspond 
to a stationary axisymmetric perfect fluid without extra symmetries and admit a 
zero pressure surface. In spite of the effort and interest put in the problem, it has 
proved difficult to obtain non-singular solutions of this kind. To our knowledge, the 
only candidate has been for a long time the Wahlquist metric jl[ . Summing to 
this difficulty, there are the ones arising from the matching with the asymptotically 
flat exterior. In fact, the derivations of the impossibility of matching Wahlquist 
with an asymptotically flat exterior [H, Isi 0] come from the analysis of the shape 
of its surface and involve approximations. This situation leads, if one is to extract 
information about these objects, to approximate solutions. 

Within the field of approximations, one has to choose between the accuracy and 
closeness to the real physical problem that numerical methods now provide and 
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the density of information and greater flexibility for theoretical work that analytic 
perturbation theory offers. Both fields have succesfully treated several important 
problems along the years in the context of stellar models and binary systems. Among 
them, the still fresh revolution in the field of coalescence of black holes or compact 
stars, a problem of fundamental importance for the modern experiments of gravi- 
tational waves detection, is actually a collaborative effort of analytic and numerical 
groups j5|-[l|. 

In ^ (CMMR from now on) some of us presented a new analytic approximation 
scheme focused on the obtention of stellar models. It is mainly a post-Minkowskian 
treatment with a secondary "slow rotation", or -to avoid confusion with the con- 
ventional use of the name- more precisely a small deformation approximation. It is 
suited to static or rigidly rotating fluids with a barotropic equation of state (EOS) 
and provides the fully matched exterior and interior spacetimes in harmonic gauge. 
The procedure is quite systematic and can be easily iterated to give better precission. 

In this paper we apply CMMR to obtain the global solution for a fluid with EOS 
e + (1 — n)p = eo, where e, p and eo are energy density, pressure and a constant with 
units of energy density, and is a free parameter. Among the relevant EOS included 
in this family are the ones corresponding to the Wahlquist fluid {n = —2) and strange 
quark matter with a simple version of the MIT Bag Model {n = A). We also show 
that the full matching procedure used in CMMR and its application to polytropes 
(lo| can be made completely general using Darmois matching conditions [ll| instead 
of Lichnerowicz ones [l2|, as well as ensuring the generality of the assumptions made 
on the embeddings of the matching surface. Later, we use the approximate metric 
to exclude Petrov type II interiors with this EOS as candidate sources of stellar 
models. This is unexpected and constrains the possible Petrov types of matchable 
rotating interiors with this EOS to I or D. It also coincides with a lack of exact 
stationary axisymmetric perfect fluid solutions of that Petrov type 3], what leads 
us to wonder if this can be the case irrespective of the EOS. This analysis shows as 
well that the stationary (non-static) constant density case can only have Petrov type 
I. Finally, we obtain the conditions for our interior to correspond to an approximate 
Wahlquist metric and recover the result that Wahlquist can not be matched with 
stationary axisymmetric asymptotically flat exteriors js], 0|- We also show that our 
interior can not be a source of Kerr metric. 



II. Interior Solution 



Our objective is to build an approximate stellar model within the context of Gen- 
eral Relativity, i.e., a global spacetime {Ai,g), where is a differentiable mani- 
fold and g a Lorentzian metric, corresponding to a stationary and axisymmetric 
spacetime containing a compact and rotating perfect fluid surrounded by asymptot- 
ically flat vacuum. In practice, this is accomplished obtaining the interior spacetime 
{V~,g~) that contains the perfect fluid and matching it through its boundary with 
the exterior vacuum solution {V~^,g^). First, we deal with the solution for the 
perfect fluid interior (y^,g^). 
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A. Source characterization 



In the interior spacetime (V~,gr~) let us call ^~ the time-like Killing vector field 
associated to the stationarity. Let also C~ be the space-like Killing vector field 
corresponding to axisymmetry, i.e., a Killing field with closed orbits and vanishing 



module on the symmetry axis|l4j]. These two vector fields commute [15 



[r,cl=o. 



(1) 



and thus two coordinates in V can be adapted to the symmetries. Let them be 

{t, (j)}, so we have 

r = dt and C = d^. (2) 

We assume the fiuid fiows on the 2-surfaces generated by the Killing vector fields so 
its 4- velocity is 

u = ip{dt + ud^) (3) 

with ip the normalization function and u the module of the rotation speed. If the 
fluid has no energy flux this is equivalent to absence of convective motion (igI. 17 . 
It also implies the veriflcation of the circularity condition [Ts], 19 1 and then the inte- 



grability of 2-surfaces everywhere orthogonal to the transitivity ones of the isometry 
group, what allows the metric to be block diagonal {Papapetrou's structure) in a 
certain set of coordinates {t, 0, r, ^}[2(|. The function u is taken here constant, so 
the fluid is rigidly rotating and thus free of expansion and shear. 

We will study a fluid with equation of state (EOS) e + (1 — n)p = cq, where e 
is energy density, p is pressure, n is a real number and eo is a constant with units 
of energy density. This linear barotropic EOS contains the ones of two signiflcant 
exact solutions: n = 1 corresponds to constant density of the Schwarzschild interior 
and n = —2 gives the EOS of the Wahlquist metric jl|. 

The assumptions of perfect fluid, circularity and rigid rotation allows us to turn 
the Euler equation VqT^ = into |2l| 



dp = {e + p) dlnip . 
It can be integrated to give for n ^ 



p 









(n-l) 


n 




£o 




n 





+ 1 



and 



(4) 

(5) 
(6) 



where iJjy: is the value of the normalization factor ijj on the surface of zero pressure. 



B. Approximations and solution 

In the remainder of this section we apply the scheme developed in Q to the 
present EOS, summarizing its main ideas. 



4 



Assuming that the gravitational field is weak enough, it makes sense to decompose 
the exact metric ga/s as 

gap = VafS + hal3 ■ (7) 

Here, rja/s is the Minkowski metric and hap, the deviation from fiat geometry, collects 
all the terms of the expansion 

ho = A/iS + A'/^S + ■ • • (8) 
hij = \h^^ + X^hij) + ■■■ (9) 
ho^=u{\h^^^ + X'h^^^ + ■■■) (10) 

where A is a dimensionless parameter related with the strength of the field. We 
choose it as A := m/ro with m := |7reoro , where tq is the coordinate radius of the 
static spherical fiuid mass, and in the term hoi we have included a global u factor 
to get a static metric when oj = 0. 

This post-Minkowskian approximation will eventually give us a solution metric in 
terms of a tensor spherical harmonics expansion. To have a cut-off for this expansion, 
we introduce a slow rotation approximation. We need a rotation-related parameter 
measuring the degree of deformation of the fluid but it must be noted that, if we 
want = t/'s to deflne almost spherical surfaces instead of cylindrical ones, we need 

~ A at least. These two requirements lead some of us to introduce in |9| 

tf-^. (11) 

the ratio between the classical centrifugal and gravitational energies (also used by 



e.g. 22|), as the slow-rotation parameter. This kind of relation is expectable if we 
want to deal with gravitationally bounded systems. It also causes that h^i terms, 
naturally possessing odd powers of Q to give the expected behaviour under reversal of 
rotation direction, have an expansion of the form hi = X^^^h^oh+^^^'^^h^^^ + 0{X''/^). 

We flx the gauge to be harmonic, and label these coordinates as {x"}. to solve 
them iteratively. Then, to get hajs we have to solve the Einstein's equations and the 
harmonic conditions system. The post-Minkowskian expansion allows us to solve 
them iteratively. The i-th. iteration gives the (9(A*) part of the metric deviation, h^^^ 
by solving the system of PDEs 



"-ka 2 ' 



where A is the flat Laplacian, Tap = T^p — \gai3T with Tap the energy-momentum 
tensor and N^p and collect quadratic and higher order terms of the Ricci ten- 
sor and the harmonic condition respectively (every rising, lowering of indices and 
operators is done with Minkowski's metric. 

The general solution of this system can be splitted into a particular solution of 
the inhomogeneous system plus the general solution of its homogeneous part, which 
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is the same at every order. Therefore, the contribution of the later to the solution 
up to order i can be written as 



1=0,2 ;^j^3 



+ J2 ^^^^r^EI + Xn% F* + 0{n^) . (13) 

1=0,2 '"O 



This expression requires several definitions and comments. First, it must noted that 
mi, ji, ai and 62 are, in general, series in A and fl (see CMMR for details). Next, 



Ti 
Zi 



Pi (cos 9) 5ijdx'® dx^ (/ > 0) , (14) 
P/(cos^) (a;*®a;^ + a;'^®a;*) (/ > 1) , 



are spherical harmonic tensors and = dt, uj^ = dr, uj^ = rd6, uj'^ = rsmOdcf) 
is the Euclidean orthonormal cobasis correspoding to the spherical-like coordinates 
{t, 0, r, ^} associated to the cartesian-like harmonic coordinates {x"}. The tensors 

e: ■■= ^ [(6 + Am - m - \{hi + 

/ 2 1 (15) 

F; :=-(/ + 1)(/ + 2)Hi - (/ + 2)Hl - -iff 

are convenient combinations of Di and these other three spherical harmonic tensor^ 



Hi 



Pi{cos9) {6ij - 3eiej)dx'^ dx^ (/ > 0) , 

P;^(cos 6) {kiCj + kjei)dx'® dx^ (^ > 1) , (16) 

P^{cos9) (kikj — 'mi'mj)dx^® dx-' (/ > 2) , 



where fc,, Cj and irii stand for Euclidean unit vectors of standard cylindrical coordi- 
nates, dp = ki dx\ dz = Ci dx^, pdcf) = nii dx^. The constant factors of the harmonic 
tensors have been chosen assigning a base dependence in VL equal to the degree of 
the Legendre polynomials involved and according to whether or not the tensors have 
spherical symmetry and thus must be present when f2 — )■ (e.g. Hq can not appear 
in that case). Also, it can be shown that oq, 02 and 62 only parametrize changes of 
coordinates. Lastly, the equatorial symmetry imposes that every summation runs 
only over even index /, except the one for Zi which runs over odd index. We cut 
the expansion at 0{VL^). 

The first order solution of (fT2|) corresponds to the constant energy density prob- 
lem already solved in |9]. It is so because from eqs. (jH) and ([gI), the first terms of 
the energy density and pressure are e ~ eo ~ 0{X) and p ~ (9(A)^. 

Defining 77 := r/ro, up to second order the interior solution is 



^ Note that the definition of E* is different from the one used in [9[. Also, := Dq has now 
spherical symmetry and the constant bo has been renamed 62- 



6 



g- = -To + Do + (mo - r/') Tq + (oq + nio - r/') Do + 
+ [mar/^Ta + (Sag + ms) r^^iDs + &2 J^o + asr/^Fg] 



2„2 



-ao-'^{n + 2) + {n + 2)S+{2 + 3n)^]To 



+ { -2ao -^{n + 2) + {n-2)S+^f + ) Do 

'5 5 21 ' 



n 



5 10 



1 



n 



- - - U^To + - - - r^^Do 



15 10 



6 



,'2 n 3 \ 2^ 

+ ( - - 802 - m2 + - - —m2n j r] D2 

I 2-rv -r 2J-^^J 21 2^ 2 y 105 

^(2mo + ao) + U2 + ^(2 - ma)) r^^ j F2 



+ X^/^Qt]^^ I ^ji - 12(ao + mo) - 3mon + QnS + 



27 15n 



7 14 



+ 



4262 A. nm2 42a2\ Sr?^ , _ 

^+Jim2+ f4-3m2-^ + 2n-^j ^ ) Zi 



4 48a2 5j3 4m2 2n 171211 > 

1 1 ri Z-x 

9 5 9 3 9 3 ' ' 



(17) 



where Ei and Fi are given in eq. fl20l ) below. They are used for later convenience 
during the matching and because the expression of g~ is shorter than using their 
interior counterparts. The constant S comes from the expansion of tl) on the surface, 
■j/^s = 1 + 'S'A + O(A^), and takes the value 



(18) 
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III. Global Solution 



A. Exterior Solution 



The approach in V"*" is entirely similar. We have to solve the system ( [T2|) for 
vacuum, i.e., with Tap = 0. Its general, homogeneous, regular at infinity solution is 



1=0,2 



(=0,2 



/=1,3 
lJ+3 



^l+l 



B. 



j,l+3 



Ei+2 + xn'4^F2 + o{n'). (i9) 



Here Tj, Di and Zi are defined as in ffT4l). and 



Er.= \l{l-l)Hi + {l-l)Hl-]^Hl (/>2) 
Fi := \l{2l -l)Di- + l)Hi-l- {Hi + fff ) 



(/>!) 



(20) 



are convenient purely spatial combinations!! The constant factors of the harmonic 
tensors have been assigned following the rules already stated for the interior solu- 
tion. Its important to note nevertheless that, since we work in spherical coordinates 
associated to harmonic Cartesian-like ones, Ml = Xn^rj+}Mi and Ji = X^n^rl+^Ji 



are the Thorne-Geroch-Hansen multipole moments |23|-|25| and An and -B„ are con- 
stants parametrizing changes within this family of coordinates. Notice also that the 
coordinates have been named as the ones used in V~ but they are not the same in 
principle. Even after the matching, the coordinates will only be C° on the surface. 
Properly speaking, {t, (p, r, 9} should be replaced by {T, $, R, 9} by now. 
The general vacuum solution is 



9^ 



-To + Do + X- 



V 



2Mo(To + Do) + 



1 



+ n'— ( 2M2{T2 + D2) 



-2Ml 



+ 



Mn 



+ 



+ ^-2 



+ I 2Mo(52 



^2 

AqMo 

^2 

2AoMo 
QB2 



B2F2 H — ttE^ 



A 
V 



9A,' 



4M2 

E2 



^2 2r]^ 



2M2 



+ 



4^4 
2M2 3^2 



Ho 



5r]^ 2r/4 
3A2M0 3A0M2 



7]^ 



7]^ 



^ Note that the definitions from Q have been modified. Now E2 has spherical symmetry. 
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+ 2Mo B: 



32Mo 



21 

2Mo 



+ 



21 

M0M2 Ao 



- 13M,) - 5^ + '-^) D, 



7 



+ 



35?7 



F2 
E4 



+ 



QB2J1 



AB2J1 ^ 9A2J1 
5 Sr/^ 



J3M0 - J1M2 + 



Zi 

3A0J3 



o(A^^]^ 



(21) 



B. Matching 

The last step to build our global spacetime M. is the matching of (y~,g~) and 
{V~^,g^). It requires first a point to point identification of a hypersurface S~ of V~ 
and another one S"*" of V'''. Later, we need to impose some conditions on g~ and 
on these surfaces in order to ensure a reasonably smooth behavior of both metrics 
considered as parts of a solution of Einstein's equations on (matching conditions). 
The above mentioned identification between hypersurfaces can be implemented by 
looking at E~ and S+ as the embeddings into V~ and of a three dimensional 
manifold E, to which we refer hereafter as the matching surface. Even though several 



matching conditions have been used in the literature [26[, it is widely accepted 
that those proposed by Darmois are the more general and suitable to the matching 
problem [llj. They impose the first and second fundamental forms of S~ and S+ 
to be equal, i.e. continuous through the matching surface. In particular, these 
conditions imply that in the neighborhood of every point of S there is a local set 



of coordinates in which the metric and their first derivatives are continuous [12 
These coordinates are called Lichnerowicz admissible coordinates and to require the 
metric to be of class on the matching surface is known as Lichnerowicz matching. 

Darmois and Lichnerowicz matching conditions are equivalent in the sense that 
the latter are the practical realization of the former in a certain set of coordinates. 
Nevertheless, if one fails to match two metrics written in some definite coordinate 
systems using Lichnerowicz conditions, it does not imply that the metrics can not 
be indeed matched. 



B.l. Matching surface 



Here, we are going to match the exterior and interior solutions given in the 
previous sections keeping all the free constants they have. First, let us start by 
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discussing how to choose the matching surface. In a general matching of spacetimes 
the surface, along with its embeddings, is actually part of the solution and can not be 
given beforehand without risk of losing generality (see jl^l and [28| for a discussion 
of this topic). In our case and in stellar model building, S~ is uniquely characterized 
in V~ as the locus of points where p = and it is the only relevant matching surface. 
The zero pressure surface r = rs(^) of our interior metric is defined implicitly by 
the equation \E'(r, 6*) = ^y; (see eq. (p)). Therefore, we may write 

S- = {t = r, = r = rs(t9), 9 = (22) 

where {r, ip, {}} are coordinates of S. 

To simplify the resolution of the matching, we are interested in using a common 
expression for S"*" and S~ in terms of the interior an exterior coordinates, that is to 
say 

S+ = {T = r, $ = y., = M^), = ^)} , (23) 

{T, $, i?, 0} being coordinates in V"*". This assumption implies no loss of generality 
if we prove that the coordinates that verify the above equations for E~ and S"*" 
belong to the class of coordinates we use to write the metrics. Then we can go along 
with the matching using the expressions for the interior and exterior metrics we got 
in the previous sections. 

Let us set up the problem. The following expression for 

S+ = {T = r, $ = i? = i?(i9), e = e(t9)}. (24) 

is well suited to the symmetries of the exterior field and how they have been im- 
plemented in the metric (!2T|) . Coordinates T and $ are adapted to the Killing 
fields, and they have been chosen to ensure that the metric tends to the flat metric 
in standard spherical coordinates at infinity. They are unique up to an additive 
constant we can set equal to zero. On the other hand, coordinates R and G are 
not completely set. Any pair of functions F{R, B) and H{R, B) leading to a set 
of cartesian-like harmonic coordinates X' = F[R, B) cos$, Y' = F{R, B) sin$ and 
Z = H{R, B) defines implicitly a couple of new coordinates by means of these two 
equations i?'cosB' = F{R,Q) and -R'cosB' = H{R,Q). Nevertheless, we must 
impose some conditions on the two functions in order to preserve the good behavior 
of the coordinates at infinity, namely 

F(i?,B) ^ /2 sin B, H{R,e) ^ RcosQ {R ^ oo) (25) 

This freedom is actually included in our metric fl2T|) by means of the constants Ao, 
A2 and B2. 

The harmonic condition requires F and H to be solutions of a second order elliptic 
system of two equations. If we add to the above mentioned boundary condition at 
infinity this other one on S"*", 

F{{Rj,{^),Qj:{^))=rj,{^)sm^, if (i?s(^), es(^)) = rs(^) cos^9 , (26) 



the solution is unique. The equation of the surface E"*" in the new coordinates (we 
drop out the primes) takes the form we want, eq. (1231 ). 
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There is another problem which should be mentioned here, even though it has no 
consequences on the equation for the matching surface itself. It has been pointed 
out that coordinates T and $ can not be naively identified with t and on the 
matching surface as we have done in the precedent paragraph [28|. Anyhow it can 
be done by making a suitable linear change, t = at' and = 0' + abt' preserving the 
regularity of the symmetry axis. 

In order to make such a change compatible with the approximate interior metric 



(eq. (llTl)). we have to assume an expansion of these two constants in powers of A as 
follows, a = 1+C(A) and b = 0{X^/^). The first and the lack of a C(A°) term in 6 just 
take into account that the starting point of our approximation is the fiat metric; the 
reason behind the semi-integer expansion of b is that under this change the angular 
velocity of the fluid reads u — b, and within the CMMR scheme this quantity must be 
of order (9(A^/^). These infinitesimal expansions actually do not keep the structure 
of the interior metric because g^^, contains a (9(A^/^) term proportional to b which 
is absent in g^^ (that starts at (9(A^/^)). It does not matter. The matching sets 
it equal to zero because the component g^^^ can not have a 0{\^^^) term unless it 
violates the asymptotic conditions on the coordinates. Moreover, the rest of the 
contribution of a and b to the metric in the new coordinates can be absorbed into 
the free constants uiq, oq and ji. 

We actually think that the choice of t and discussed above has already been 
considered in the CMMR scheme from the very beginning. Dropping out a (9(A^/^) in 
g^^ at the linear level of the approximation is a way of choosing the inner coordinates 
t and 0, and also the angular velocity w in a sense, since we asumed g^^p to be 



proportional to u in eq. (jlOl ). 

These kind of assumptions are meaningful in a perturbation scheme but they 
can not easily be implemented in a exact matching problem. The continuity of the 



Killing fields on the surface matching used in [28| is a smart reasonable assumption to 
make an exhaustive use of the symmetries of the problem. The arguments sketched 
above show how this point of view has not been forgotten in our scheme. Therefore, 
we can argue that our approach is in accordance with it. 

Lastly, let us simplify a little bit more the matching process by introducing an 
ansatz for the explicit equation of the surface. Being this an axisymmetric problem, 
the zero-pressure surface can be expanded as a power series in Legendre polynomials, 
and the reflection symmetry allows us to rule out every odd degree term in the 
expansion. Coherently with the assignment of dependence on Q we made for the 
approximate solutions of the homogeneous linear Einstein equations, we have then 

MO) = ro [1 + ^]VP2(cos 6)] + Oin^) (27) 

where tq is the coordinate radius of the fluid at rest and a a constant to be deter- 
mined while matching. 

B.2. Darmois matching 

First we impose Darmois conditions on the surface fl27|) to match {V^,g^} and 
{V~^,g^} up to 0{X'^,Q^) in the metrics. They are satisfied when the constants 
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associated to the multipole moments are 
'3ao n 14 2 



Mn 



1 + A 



1 



Ms = -- + A 
Ji 



2 

2 



+ A 



ao 



2n 16 

35 ^ y 



5ao 
"6~ 



3n 176\ 

35 ^ Tosy 



1 



^ + A( + 



Oq 11^ 



2 441 735 J 
while their interior counterparts take the values 
3n 9 



(28) 
(29) 
(30) 
(31) 



mo 
ji 



3 + A 



3ao H ^ h 1 

"42 V 2 



+ C(A^fi^ 



-l + A(-ao + 262-^-|~)+0(A^fi^), 



2^]2 



+ A 



49 ^2 

4ao H ^ h 1^ 

" 2 5 



1^ + 26,-^ 
3 14 



289 
105 



4ao ^ ^ ^ 1262 
+ 2a2 + 



7 25 
and finally, the exterior gauge constants are 



^_^)+o(A^^]^^ 

49 245 ' ^ 



(32) 
(33) 

% 
(34) 

(35) 



^0 = ao + A 



3a2 



^2 = -y + ^2 + A 



So 



f +62+A 



3ao — 



237 22^]2 



35 



35 



7a2 



3a2 



+ 0200 + 

25ao 



+ C(A^^] 
139ao 



70 



— 3a5 



362 - - 
2 



383 
"90" 

i2 



o(A^^]^ 



C(A^fi'). 



7 8 14 

The matching surface on which the matching conditions hold is 



3a2 , 5?2 10 

— ^ + &2 H \ 

2 21 21 



c(A^^]^). 



(36) 
(37) 
(38) 

(39) 



These results require some comments. Here, the only free parameters are ao, 02 and 
62, (which, as already mentioned, parametrize changes of harmonic coordinates in 
V~), To (which depends on the size of the source hut also on the coordinates), eo, 
a; (which are part of the definitions of A and VL) and the EOS parameter n. Then, 
for a fixed set of source parameters s := {n, eo, uj} and roH the interior metric that 



^ Note that has been exckided from s because of its coordinate dependence. It contains infor- 
mation necessary to fully characterise the source, though. 
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can be matched is unique up to changes of coordinates. Since the asymptotically 
flat exterior of a certain source spacetime in rotation is unique 28| then, given a 



set of source parameters s and tq, there is only one global spacetime possible, i.e., 
only a couple of metrics {g~ , g^) among the families g~ and gr"*" give spacetimes 
that can be matched. Nevertheless, this could seem contradictory with the apparent 
fact that the value of mass and angular multipole moments depend on the value of 
ao as eqs. (28) to (jsil ). This apparent dependence happens because these constants 
are not the only gauge dependent quantities in the expressions. There is coordinate 
dependence hidden in A, which, unlike f2, depends on tq. 

This problem can be solved finding their expression in terms of physical constants 
because they are gauge invariant. A convenient way of doing it is using the mass 
monopole moment Mq to redefine A. In CMMR, A is defined as A := ^eo'^o ^^^^ 



Mn 



— lim 



]9tt,r 



Xr-oMo 



(40) 



and Mq = 1 + 0{X) refiecting the fact that A and tq were chosen to reproduce 
the Newtonian mass of the source. Now, if we define {A', Tq} to give the GR mass 
monopole moment Mq and keep the same relation between them, i.e. 

4:71 ,2 



A': = 



AVn 



Mo, 



(41) 
(42) 



inserting eq. fl4ll ) into eq. f l42l ) and using the expression of Mq 
XtqMq, we obtain that 



to find Mn 



A' = A ( 1 + A <( ao + ^ 



- H \ Q^(4 - n) 

5 5 15 ^ ' 

n 14 2 

5+T + T5^(^--) 



^o(A^^]^), 
+ o{x\n^) 



(43) 
(44) 



Using these changes, Mq = 1 + C(A'^) as required and the dependence on the gauge 
constant ao of the multipole moments M2, Ji and J3 dissapears. In this way, the set 
of parameters to completely specify the interior would become eo, n, u and Mq. 

This procedure can be followed as well using the central pressure Pc instead of 
Mq to characterise the interior as is sometimes done in astrophysics. In this case, 
the new approximation parameter A is 



A = X<1 + X 



with 



3n 17 ^9/3 n 



A := p,- ( 1 + -n^ ) + 0{Q 



+ c(A^^]^ 



(45) 



(46) 



Keeping a relation of the form of eq. (1411 ). the associated radial coordinate changes 
as 

ao 'in 17 ^r, / 3 n 

— H 1 \ 

2 10 10 VlO 12 



r-A = ro < 1 + A 



0(A^fi^). 



(47) 



13 



Again, as one expects, with these changes the expressions for the multipole moments 
become manifestly coordinate independent. 

Thus, with these last results, the exterior ( 1211) and the interior (1171) with their 
constants taking the values eqs. ( 1281 ) to (1381 ) give the most general approximate 
family of global asymptotically flat solutions for the kind of source studied, each of 
its members characterized only by the values of {n, eo, oj} and Tq or r\ that, with 
the first three fixed, depend only on Mq and pc, respectively. Additionally, they also 
point out the behaviour one intuitively expects as a generalization of the theorem in 
[29|, i.e., that for a stationary axisymmetric singularity free compact rotating perfect 
fluid, its asymptotically fiat exterior is unique once the EOS, central pressure and 
rotation speed are fixed. 



B.3. Lichnerowicz matching 



Now we impose Lichnerowicz conditions . The 0(A^, fi^) metric is then matched 
in the global set of coordinates when its multipole moments and {mi, ji} constants 
take the values 



Mo 
M2 
Jl 
h 
mo 
m2 
jl 



1 + A 
1 



+ A 



n 14 ^0/8 2n 

-H hfi 

5 5 1 

n 37 
14 ^ 35 



15 15 
+ 0(A^^]' 



+ C(A2,fi^ 



2 fi2 
5 + T + ^ 



2^ 16 , ^2 



4- 



35 7 
lln 496 



/176 


3nY 


V105 " 


357 



441 735 



3 + A 



3n 9 r^o n 

\ hfi l 

4 2 V 2 



0(A2,fi2), 

fO(A2,fi^), 



2r]2 
2 + -- + A 



^ 49 ^2 

- H h fi^ 

2 5 



289 
105 



5n\ 



+ o(A^^]^), 



7 



+ A 



"in 
'49 



326 A 
2457 



+ 0(A^^]2 



and the coordinate-parametrizing constants are 



4A 



63 

B2 = o{\\n''), 



A 7 + 



2^]^ 



+ c(A^^]^ 



86A 



02 = 

.2 = o(A^^] 



105 



+c(A^fi' 



(48) 
(49) 
(50) 
(51) 
(52) 
(53) 
(54) 
(55) 



(56) 



The value of these parameters is unique for each set of parameters {n, eo, tq] 
and verify the relations obtained with Darmois conditions as is to be expected. 
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Nevertheless, it proves the existence of a harmonic and asymptotically Cartesian 
global system of coordinates up to the approximation order considered. 

In spite of the attention drew to Darmois matching before, it is important to 
remark that the real focus of the approximation scheme is the obtention of totally 
matched spacetimes. In fact, only Lichnerowicz conditions were used in 10 1, al- 
though the generality of the results in them can be verified with the same techniques 
used here. Besides, one needs the matched spacetime for many practical pur poses, 
as for example to compare with numerical results for stellar models jsol . 3ll | built 
using global coordinates. 



IV. Petrov classification 

We will now analyse the possible Petrov types of the unmatched interior metric. 
They are given by the possible Jordan canonical forms of the Q"^ matrix, defined 
from the Weyl tensor C"^^^ as 

Q"^ := (C%;. + i*C%M)MV 

where * denotes the left Hodge dual and is a unit timelike vector. We use m* = 
for simplicity. Due to the symmetries of the Weyl tensor, Q^^m" = 0. This implies 
that for the eigenvalue problem 

Q"^t;" = ev'^ (57) 

v"Ua = and thus, for a metric with Papapetrou's structure, ^ is completely 
classified studying 

Q=\Q\Q\ . (58) 

V Q\j 

The matrix Q is always trace-less and in an orthonormal cobasis also symmetric. 
For computation convenience we work now in spherical-like coordinates and thus 
6 7^ r general. The eigenvalues e\ and of the r — Q subspace are degenerate 
iff the discriminant of the roots of its characteristic equation is zero , i.e., 

[Q\-Q\)' ^^Q\Q\ = ^ (59) 
which using the traceless property gives 

Q\Q\- Q\Q' r-\[Q\f = (60) 

The condition for either e\ or £2 to be degenerate with £3 = Q'p ^ is 

Q\Q\-Q\Q\^m\f = ^. (61) 
that has the same expression as eq. fjool ) switching the numerical factor. The struc- 



ture of eq. ( 15 8l ) allows the following possibilities. First, no degeneracy at all. This is 
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the general case and corresponds to Petrov type I. The second, degeneracy in only 
two eigenvalues. It can come from eq. (j60l ). in which case it can lead to types D 



and II, or from eq. ( l6ll ) that can only give rise to type D. Last, three degenerate 
eigenvalues, what happens if both eqs. B) and Q are fulfilled. This, with the 
traceless property rules out type III and hence the only possibilities are types N and 
O. 

We get the following conditions on the metric parameters (in which the trivial 
case A = is not considered) extracting from (160|) and fl6T!) the relevant information 
up to this order of approximation. 

We will first analyse the general rotating case separately from the static case for 



clarity. The relevant condition eq. (1591 ) for ei =62, even considering only its terms 
up to C(A3,fi3) 

-X^/^Q^^im2Pi + X^Q^^ [5m2{n - 1) (4P2 - 1) - 108] + O{X'^\Q^) = (62) 

can not be satisfied everywhere unless = and therefore we can not have ei = 62 
degeneracy out of the static case. The conditions for the other possible degeneracies 



El = £3 or 62 = 63 are given by eq. fl6ll). that yields 

-X'n'^{P2-l) [5m2(n - 1) + 18] 

+ iX'/^n'-^ (Pi - P3) [2 (35j3 + 3) (n - 1) + m2(23 - 14n)] + 0(A^ n") = 0. 
175rU 

(63) 

They are satisfied in the static case and, if 7^ 0, when the constants of the metric 
verify 

18 , . 3(n + 8)(8-5n) 

m2 = 37^ r and J3 = — 7^777 tt^ — , 64 

5(1 — n) 175(n — ly 

which can never be satisfied in the constant energy density case. 

Concerning the static case, equation ( 15^ gives conditions different from eq. (}6i 
They are satisfied only when n = 1, while the other degeneracy possibility condition 
0611) is always verified. This can be seen straightaway from the form the Q matrix 
takes 

. 2 /2(n-l) \ 
Q = -^aM 1-n ] (65) 
^^0 \ 1-n 



as well as the fact that the only possible Petrov types are D(n 7^ 1) and O (n = 1) as 
must be the case for a spherically symmetric spacetime |32l . p. 228]. The condition 
for type O is expectable since an asymptotically flat static solution with our sym- 
metries must have constant energy density [33^]. From the equations fl59|) and fl6T]l 
alone one can not say whether this behaviour will endure further approximations 
because extra conditions could impose a more general Petrov type, but these theo- 
rems on spherical symmetry and constant energy density give a strong hint about 
its endurance. 
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Therefore, collecting the results together from the rotating and static cases, we 
conclude that 

• An invariant subspace of dimension 3 for the eigenvalue problem, which would 
lead to Petrov type III, is ruled out by the structure of Q associated to the 
Papapetrou structure of the metric. From our perturbation theory results, 
we see that the only option for a bidimensional invariant subspace appears in 
the static limit, where its existence is forbidden by the fact that the spherical 
symmetry associated imposes types D or O. Accordingly, the Q matrix of our 
interior spacetime is always diagonalizable. 

• In general, the Petrov type is I. Out of the static case, it will only be type D 



when eq. (l64l ) are satisfied provided n ^ 1. In the static case, it will always be 
type D unless n = 1, in which case the Petrov type is O. Then, the constant 
energy density case can only be type I(f2 7^ 0) or 0(f2 = 0). 

It could be argued that our approximate results do not necessarily hold for 
exact solutions. Nevertheless, as long as an exact solution for the source we 
work with exists, a series development of it following the CMMR approach 
must lead to our results. Because of this, while any property compatible 
with a truncated series development is not necessarily a property of the exact 
solution, a behaviour ruled out already in the truncated solution can not be a 
property of the hypothetical exact solution. 



V. Some implications 



It has long been suspected that the Kerr metric can not represent the exterior of 
any stellar model. We can easily check that it is indeed the case here using the kind 
of analysis that appears in j9|. 



The first three Kerr multipole moments are [25 



Kerr 



Kerr 



ma. 



Kerr 



—ma 



(66) 



If our two first multipole moments were equal to the Kerr ones, M^*^^^ should have 
the expression 



Kerr 



r 



m = XtqMq 
ma = X^/^Vtrl-h 



-ma 



XroMo 



i.e., the first A-order component of M2 should vanish. This in in contradiction with 
eq. (j49l ) and hence neither our interior nor any exact metric of which it could be an 
approximation can be a source of Kerr. 

Now we focus on the Wahlquist family of metrics [l[ . That is to say to a stationary 
axisymmetric rigidly rotating perfect fluid with EOS e + 3p = const, and Petrov type 
D [SJ, |35|. In the non-static case, the conditions for Petrov type D for our interior 
family are given by eq. (64)- Then, in the case when n = —2 and eq. fj63 ) is satisfied, 
g~ is of type D and our interior is an approximation to the Wahlquist family. 
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It must be noted nevertheless that, despite the fact that Wahlquist's family is a 
subcase of our general interior solution (as must be expected), the values m2 and p 
imposed by the matching with the asymptotically fiat exterior do not satisfy eq. (l64l). 
Then, we recover and give an independent derivation of the known result that, within 
perturbation theory, Wahlquist's family can not correspond to an isolated source 



A last comment. If one looks for stationary axisymmetric perfect fluid solutions 
with a static limit as candidates to represent the interior of a stellar model, consi- 
dering the Penrose chart of how a certain Petrov type can lead to another through 
degeneration, one should then only take into account those metrics whose Petrov 



type can lead to the types D and O corresponding to spherical symmetry [13 



A type N, rigidly rotating perfect fluid with barotropic EOS and e + p ^ can 



not be axisymmetric |36[, therefore we must discard types III and N, but all the 
rest should, to the best of our knowledge, be considered. It seems reasonable that 
a type II exact metric with the properties we demand can be approximated by 
our solutions. Nevertheless, the Petrov type II is not included among the possible 
types of our general interior metric and hence, we claim that there is no stationary 
axisymmetric rigidly rotating perfect fluid metric witli EOS e + (1 — n)p = const, of 
type II possessing a static limit and a surface of zero pressure. This is in accordance 
with the weird fact that, even dropping the demand of zero pressure surface, it has 
not been found any type II exact interior metric suitable to be part of a stellar 
model, while the harder field of type I solutions is not empty [13(| . 
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A. 0(A^, Q^) metric components after Lichnerowciz matching 



Here we give both metrics written in the orthonormal cobasis associated to 
It, r, 6, d)}. These result from the substitution of eqs. (jssl ) and (56) in eqs. (17) 
and (12 ll). The exterior components are 
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The results for the interior metric are 
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